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We study systems of particles on a line which have a maximum, 
are locally finite and evolve with independent increments. "Quasi- 
stationary states" are defined as probability measures, on the a- 
algebra generated by the gap variables, for which joint distribution 
of gaps between particles is invariant under the time evolution. Ex- 
amples are provided by Poisson processes with densities of the form 
pidx) = e~ ax sdx, with s > 0, and linear superpositions of such mea- 
sures. We show that, conversely, any quasi-stationary state for the in- 
dependent dynamics, with an exponentially bounded integrated den- 
sity of particles, corresponds to a superposition of Poisson processes 
with densities p(dx) = e~ ax sdx with s > 0, restricted to the relevant 
er-algebra. Among the systems for which this question is of some rele- 
vance are spin-glass models of statistical mechanics, where the point 
process represents the collection of the free energies of distinct "pure 
states," the time evolution corresponds to the addition of a spin vari- 
able and the Poisson measures described above correspond to the 
so-called REM states. 

1. Introduction. Competitions involving large numbers of contestants 
are an object of interest in various fields. One could list here the energy 
levels of complex systems and the free energies of competing extremal states 
of spin-glass models [10] and include a broad range of other examples. We 
are particularly interested in dynamical situations where the competition 
continues in "time," though time may be interpreted loosely. For example, 
in the motivating example of spin-glass models [10], a point process on the 
line represents the collection of the free energies of distinct "pure states" of 
a system of many spin variables, and the "time evolution" corresponds to 
the incorporation in the system of yet another spin variable. 
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Influenced by the terminology of statistical mechanics, we use here the 
term state to mean a probability measure on the relevant cr-algebra of subsets 
of the space of the point process configurations. For much of the discussion 
which follows, the relevance would be limited to the information concerning 
only the relative positions of the points, relative to the one which leads at 
the given instant. 

As in the pictures seen in marathon races, often the point process de- 
scribing the relative positions appears to be time invariant. We refer to such 
states as quasi-stationary. 

In this paper we characterize the quasi-stationary states for the class of 
systems in which the evolution occurs by independent identically distributed 
increments of the individual contestants. The main result is that any such 
state, of a point process with locally finite configurations with more than 
one point and exponentially bounded density, corresponds to a linear super- 
position of Poisson processes with densities of the form 



with s > 0. This may be rephrased by saying, in the terminology coined by 
Ruelle [11] (who invokes the work of Derrida [6]), that all quasi-stationary 
states correspond to superpositions of the random energy model (REM) 
states. 

Remark. Our main result may have a familiar ring to it, since the 
above distributions are known to describe the "Type-I" case of the extremal 
statistics [8]. 

Remark. It would be of interest to see an extension of the classification 
of the quasi-stationarity to a broader class of dynamics where the evolution 
may exhibit correlations. One may note that the REM states have an ex- 
tension, based on a hierarchical construction, to the family of the GREM 
states [2, 11], which exhibit quasi-stationarity under a broad class of cor- 
related dynamics. Is that structure singled out in some way by its broader 
quasi-stationary properties? 

In the following section we introduce the concepts more explicitly. We refer 
to the system as the Indy-500 model, ignoring the fact that for a number 
of obvious reasons this is not a faithful description of the dynamics in that 
well-known car race. 

2. The Indy-500 model. The configuration space of the Indy-500 model 
is the space fi of infinite configurations of points on the line, which are 
locally finite and have a maximum (in the order of R). Its elements wGfi 
can also be described as sequences, uj = {xn}ra=i,2,..., with 

(2.1) x\ > X2 > ■ • ■ and x n — > — oo. 




p(dx) = e sx sdx 
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(Variables written as x n should by default be understood to be ordered.) In 
the time evolution considered here the points evolve by independent incre- 
ments. 

As is generally the case with stochastic evolutions, the dynamics can be 
presented in two ways: as a stochastic map, in which the configuration lo G Vt 
is changed in a random way — through the independent increments, or as a 
reversible transformation taking place in a larger space, which encompasses 
the full information about both the future and the past dynamics. Our 
terminology is based on the former view; however, the second perspective 
provides a useful basis for the intuition guiding the analysis. 

Thus, the time evolution is given by a stochastic map determined by the 
collection of random variables 77 = {/in}n=i,2,... : 

(2.2) T n : {x n } h-> {x n } with x n = x Un + hn n , 

where h n are independent random variables with a common probability 
distribution g{dh) on R, and II is a permutation of N, which depends on 
both u) and 77, aimed at recovering the monotonicity for x n . In other words, 
IT = IL(u,rj) is a relabeling of the moving particles according to the new 
order. 

For a given probability measure fi(duj) on O, we denote by T/j the corre- 
sponding probability distribution of the one-step evolved configuration {x a }. 
To be explicit: the average over T/j corresponds to averaging over both \i 
and 77. 

One needs to pay some attention to the <r-algebras on which the measures 
fi and Tfi are to be defined. Since we are interested in the classification 
of states which are only quasi-stationary, we allow those to correspond to 
probability measures defined on a smaller cr-algebra than the one usually 
used for point processes on a line. (Such a change makes the result only 
stronger.) 

The standard <r-algebra, which is natural for the state space of particle 
configurations, is generated by the occupation numbers of finite intervals 
(see, e.g., [4]). Let us denote it by B. Measurable functions include all ip : $7 — > 
R of the form 

(2.3) 1>f(u):=J2f(xn) 

n 

with bounded measurable functions / : R — ► R, of compact support. How- 
ever, in this work we are interested in probability measures on the smaller 
ex-algebra B generated by functions which are invariant, ip(SbUj) = ip(u>), un- 
der the uniform shifts 

(2.4) S b : {x n } 1 ^ {x n } with x n = x n + b. 
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Functions which are measurable with respect to B depend on the config- 
uration only through the sequence of the distances of the particles from the 
leading one: 

(2.5) u n = xi- x n . 

Thus, a probability measure \i on (0, B) is uniquely determined by the "mod- 
ified probability generating functional" (MPGFL) 

(2.6) G At (/)=E^exp|-^/(x 1 -x n ) 

with /(•) ranging over smooth positive functions of compact support. [The 
regular "probability generating functional" is defined without x\ in (2.6).] 

One can now formulate a number of distinct "steady-state" conditions, 
where the term state refers to a probability measure on a suitable a-algebra, 
which is not always the same. 




Definition. A stationary state is a probability measure n(duj) on (f2,i3) 
which is invariant under the stochastic map T, that is, T[i = fi, or more 
explicitly, 

(2.7) E^(Tu)) =E^(u)) 

for any ,8-measurable ip, where the expectation functional E /t includes an 
average over both u (distributed by //) and T [determined through {h n }, as 
in (2.2)]. 

A steady state is a probability measure fi(duj) on (fl,B) for which there 
is a nonrandom V (= the "front velocity") such that T/i = Sy^, that is, 

(2.8) E^(T{x n }))=K^({x n + V})) 

for all i3-measurable functions tp. 

A quasi- stationary state is a probability measure fj,(du) on the a-algebra 
B (sub-u-algebra of B) such that (2.7) restricted to shift-invariant functions 
ip holds, that is, for which 

(2.9) E M (y>(K})) = E M (v>(K})) 

with {u n } the gaps defined by (2.5), and {u n } the gaps for the configuration 
Cj = Tuj. 



For an alternative characterization of quasi-stationary measures, in terms 
which are more standard for point processes, let us note that each configu- 
ration is shift-equivalent to a unique element of the set 

(2.io) n o = {{x n }\ Xl = 0}. 
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The "normalizing shift" S:ujh^ S— X1 ( U )U induces a measurable map from 
(O, B) to (Q , B) C (f2, B), and thus also a map (for which we keep the symbol 
S) which associates to each probability measure ji on (£l,B) a probability 
measure Sp on (Q,,B), supported on . The measure p is quasi-stationary 
if and only if the corresponding measure Sp is invariant under ST — the time 
evolution followed by the normalizing shift. 

Stationarity is a special case of the steady state, and the latter reduces 
to it when viewed from a frame moving at a fixed speed. Quasi-stationarity 
is the less demanding property of the three mentioned above, and is the 
condition of interest if one follows only the relative positions. 

Through a combination of the results in [9] and [3] one may conclude that 
any steady state of the Indy-500 model, whose jump distribution satisfies the 
nonlattice condition (meaning that its support is not contained in any set 
of the form a + VL C R), is a Poisson process with a density of the form 
p(dx) = se~ sx dx. These are the REM states which are discussed in the next 
section. 

Our main result is that for the infinite systems discussed here quasi- 
stationary probability measures can only be linear superpositions (as prob- 
ability measures) of the above steady states restricted to B. 

Remark. The restriction, in the above statement, to infinite number of 
particles excludes the trivial example of a quasi-stationary state which is 
not the projection of any steady state, which is provided by a single point 
moving on the line by independent increments. In this case the state looks 
stationary from the perspective of the "leader": there is always just one 
point, at the origin. There is, however, no steady velocity V such that (2.8) 
holds. 

Remark. Linear superpositions (of measures on the suitable a-algebras) 
preserve the property of quasi-stationarity though not that of steady state — 
due, in the latter case, to the possible variation in the front velocities. 

3. The REM states. We recall that for a probability measure p(dx) on 
R, a Poisson process with the density p is a probability measure on (O, B) 
for which the occupation numbers for disjoint sets Acl form independent 
random variables, £(A;lu) = £(A), with the Poisson distributions 

(3.1) Prob(£U) = k) = ^^ e - p{A) and mean E(£(A)) = p(A). 

We denote by p SyZ {du) the Poisson process with density p StZ (dx) = se ~ s ( x - z ) d x 
on R. 

The special role of the above states in the present context is already seen 
in the following statement, which is based on known results. 
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Proposition 3.1 (Based on [[3, 9]-[ll]]). For any nonlattice single-step 
probability distribution g{dx), the collection of the steady states correspond- 
ing to the evolution by i.i.d. increments {h n } with the distribution g{dh), as 
described by (2.2), consists exactly of the probability measures p s ^ z (duj) [on 
(Q, B)], with s > 0, z G R. For each of these states, the corresponding front 
velocity V is the solution of 

(3.2) e sV = J e sx g{dx). 

Furthermore, with respect to p S)Z (du)), the past increments also form an i.i.d. 
sequence, however with a modified distribution: conditioned on {x n }, the 
variables {hu n } form a sequence of i.i.d. variables with the probability dis- 
tribution 

(3.3) Sm - 



ke s yg(dy)' 

Thus for these steady states the distribution of the increments changes 
depending on whether one looks forward or backward in time (!). In other 
words, the permutation n n (a>) transforms the sequence of i.i.d. variables 
{h n } into an i.i.d. sequence ({^n n }) with a different distribution. (Of course 
this is possible only in infinite systems.) 

Proof of Proposition 3.1. The evolution by independent increments 
is well known, and easily seen, to take a Poisson point process into another 
such process with the density modified through convolution (p i— ► p * g). 
Therefore, just the steady-state property of the states p, s>z is an elementary 
consequence of the behavior of the exponential density under convolutions. 
However, for the more complete statement made above it is useful to ap- 
preciate the following observation, concerning two possible ways of viewing 
the collection of variables u) = {(x n , h n )}. The following are equivalent con- 
structions of a point process in R x R: 

(i) A collection of points {x n } is generated via a Poisson process on R, 
with the density p(dx), and then to each point is attached, by its order on R, 
a random variable {h n }, taken from an i.i.d. sequence with the distribution 
g(dh). 

(ii) The configuration is generated directly as a Poisson point process in 
R x R, with the two-dimensional density p(dx)g(dh). 

The transition of the perspective from (ii) to (i) requires only the second 
factor in the product measure on R x R to be normalized f K g(dh) = 1. 

Now, the map (x, h) ^> (x + h, h) = (x, h) takes the Poisson process de- 
scribing u> into another Poisson process on R x R, which yields the joint 
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distribution of the "new" positions paired with the steps "just taken." In 
case of p StZ (dx) x g(dh), the density of the new process is: se~ sx dxg(dh) = 
se -s(x-h) dxg[dh). This can also be written as a product [J e sy g(dy)]se~ sx dx x 
-4 — 9 ^ dh \ , where now the second factor is properly normalized. By the pre- 
vious observation it immediately follows that: 

(i) The positions after the jump {x n } are distributed as a Poisson process 
on R with the modified density £(cfcc) = [/ e sy g(dy)]se~ sx dx = se~ s ( x ~ v ^ dx, 
that is, {x n } have the same distribution as {x n } + V with V satisfying (3.2). 

(ii) When conditioned on the configuration {x n }, the jumps just taken 
are generated by an independent process on R with the probability density 
given by (3.3), as claimed. 

For the converse statement, that is, to prove that all steady states are 
of the REM type, one may first note that if p{duo) is a steady state for 
the dynamics corresponding to g(dx), with the front velocity V, then p is 
stationary under the dynamics corresponding to a shifted single-step dis- 
tribution: gy(dx) = g(d(x — V)). The classification of stationary states, and 
hence also steady states, is found in [9] , where it is implied that any station- 
ary measure is a superposition of Poisson processes whose random density 
solves the equation p = p* g. As established in [3], for nonlattice g(dh) the 
only solutions in the space of positive measures are measures of the form 
£,(dh) = [Ae~ sh + B] dh. The condition that the typical configurations be 
bounded on the positive side imply that s > and B = 0. □ 

Having introduced the REM states, we are ready to formulate the main 
result. 

4. Classification of quasi-stationary states. 

Definition 4.1. A probability measure p on f2 is g-regular if for al- 
most every Tuj = {uj,{h n } n }, with respect to p(du>)Yl n£ ig(dh n ), the point 
configuration {x n + h n } n is locally finite, with a finite maximum. 

The ^-regularity of p means that with probability 1 the configuration 
obtained through the independent increments has a maximum and can be 
ordered. This is a preliminary requirement for the possible quasi-stationarity 
of p. It is easy to see that a sufficient condition for ^-regularity is met in the 
situation discussed next. The general sufficient condition is the finiteness, 
for all x e R, of 



(4.1) E(card{n : x n + h n > x}) = E M Prob g (/i n >x — x n ) 
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In the following, to simplify the exposition and to avoid confusion we will 
always assume that at t = the rightmost particle in the configuration is at 
x\ = (we can do this without loss of generality); we will denote by x n the 
positions of the particles at t = and by y n the positions of the particles at 

t = T. 

Following is our main result: 

Theorem 4.2. Let g be a probability measure with a density on R and 
let fj, be a probability measure on (Q,B), satisfying 

(4.2) J e sx g(x) dx < oo VsGK 

and 

E„({{j of particles within 

(4.3) 

distance y of the leading particle}) < Ae Xy \/y > 

for some A > and A < oo. If [i is quasi- stationary with respect to the 
dynamics corresponding to independent increments with the distribution g, 
then it is supported on Poisson processes with densities se~ sx dx, s > 0. 

The meaning of the theorem is that the probability space £1 can be split 
into pieces and the process on each piece of is a Poisson process with a 
density se~ sx dx for a particular s. 

In the proof we shall use the fact that point processes are uniquely deter- 
mined by their probability generating functionals (as discussed in [4]). Our 
derivation of Theorem 4.2 proceeds along the following steps. 

1. First we note that any quasi-stationary state can be presented as the 
result of evolution of arbitrary duration (r) which starts from a random 
initial configuration, distributed by the given quasi-stationary state, and 
evolves through independent increments. 

2. Analyzing the above dynamics, we show that for large r the resulting 
distribution is asymptotic to Poisson processes with the corresponding 
(evolving) densities. Thus, it is shown that the quasi-stationary measure 
fj, can be presented as the limit of a superposition of random Poisson 
processes, where the randomness is in both the Poisson measure and the 
resulting particle configuration (Theorem 5.1). 

3. Applying a result from the theory of large deviations (Theorem A.l), 
and some compactness bounds which are derived from quasi-stationarity, 
we show that the quasi-stationary measure admits a representation as 
a random Poisson process, whose Poisson densities (F) are the Laplace 
transforms of (random) positive measures (Theorem 6.1). Furthermore, 
in this integral representation of (J,, F may be replaced by its convolution 
with g, followed by a normalizing shift. 
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4. For the last step we develop some monotonicity tools (Section 7), for 
which the underlying fact is that under the convolution dynamics the 
Laplace measures increase their relative concentration on the higher val- 
ues of the Laplace parameter (Theorem 7.3). This corresponds to the 
statement that unless the function F is a strict exponential, under the 
convolution dynamics the function F becomes steeper, and the distribu- 
tion of the gaps is shifted down. Using a strict monotonicity argument, 
we show that quasi-stationarity requires the measures in the above super- 
position to be supported on pure exponential functions (or, alternatively 
stated, functions whose Laplace measure is concentrated at a point). 

The final implication is that the quasi-stationary measure is a superposi- 
tion of REM measures, as asserted in Theorem 4.2. 

Let us remark that Section 7 may be of independent interest. It is noted 
there that within the space of decreasing functions which are the Laplace 
transforms of positive measures on [0, oo), convolution with a probability 
measure makes any function steeper, in the sense presented below, except 
for the pure exponentials on which the effect of such a convolution is only a 
shift. 

5. Representation of /x as a random Poisson process. 

5.1. "Poissonization" — the statement. Let J- be the space of monotone 
decreasing, continuous functions F:R — ► [0, oo], with F(x) — > for x — ► oo 
and F(x) — ► oo for x — > — oo. We regard a function F £ T as normalized if 



and denote by M the normalizing shift: M:F(-) i— > F(- + zf), with zp = 
sup{z 6 R : F{z) > 1}. 

For each F 6 the Poisson process on R which corresponds to the mea- 
sure (-)dF will almost surely exhibit a configuration which can be ranked 
in the decreasing order of R. The probability that there is no particle above 
x € R is exp(— F(x)). Conditioned on the location of the leading particle 
(x), the rest are distributed by a Poisson process on (— oo,x] with the den- 
sity d(—F). Thus, the MPGFL [defined in (2.6)] of the Poisson process with 
density F, which we shall denote by Gf{/), is given by 



(5.1) 



F(0) = 1 




f(*-y))d(-F(y))y 



Let us note that 



(5.3) 



G F (f) = G MF (f), 



since the probability distribution of the gaps is not affected by uniform shifts. 
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For the purpose of the following theorem let S T be a random variable with 
the probability distribution P(S T > y) = J I[J2yj > y]g(yi) ■ ■ ■ g(y T ) dy! ■ ■ ■ dy T . 
We associate with each configuration u>, and r G N, the function 

(5.4) 

•E ) 3 

m 

and denote by z UT the position at which 

(5.5) F^ T (z^ T ) = 1. 

One may note that F w - T (x) is the expected number of particles on [x,oo) 
for the configuration which will be obtained from oj after r steps of evolution 
with independent increments. If the support of g(y) is not bounded, one may 
easily find configurations for which F w - T (-) diverges. However, if the measure 
[i is ^-regular, then a.s. F U . T (-) < oo. Furthermore, we shall see that if fi is 
quasi-stationary, then the position of the front after r steps can be predicted 
to be in the vicinity of z w - T — up to a fluctuation whose distribution remains 
stochastically bounded (i.e., forms a "tight" sequence) as oo. 

The main result of this section is: 

Theorem 5.1. Let [i be a g-regular quasi-stationary measure, for the 
independent evolution by steps with some common probability distribution 
which has a density g(u). Then for every positive function f of compact 
support in 1R ; 

G M (/)= lim / n{dw)G MFu , iT {f) 

(5.6) 

= lim 

where G^(f) is the modified probability generating functional defined in 

(2.6) . 

This statement implies that the measure /i is, in the "weak sense," a 
limit of random Poisson processes, of measures corresponding to the random 
functions NF^ T {-) whose probability distribution is induced from /i through 
their dependence on lj. 

Let us note that this result is related to — but not covered by — the known 
statement that any limit of a sequence of point processes which is derived 
through successive random independent increments is a mixed Poisson pro- 
cess (e.g., [4], Theorem 9.4.2). Unlike in that case, the time evolution consid- 
ered here incorporates shifts according to the position of the leading point 
(and the limiting process is not stationary under translations). 

The rest of this section is devoted to the proof of this assertion, for which 
we need some preparatory estimates. 

First let us make the following observation: 



/j,(du;)G g *tf F „. T {f) 
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LEMMA 5.2. Any quasi- stationary measure is supported on configura- 
tions with either exactly one particle, or infinitely many. 

Proof. The statement is a simple consequence of the spreading of the 
probability distribution of the sum of independent increments, that is, of 
the variable S T . For example, one may consider the function 

(5-7) Y^(y)= f ,({y 1 -y 2 >y}). 

(2) 

By the dominated-convergence theorem, Y} (y) — > 0. However, for any 

y->oo 

finite number of particles, the probability that after r steps the smallest gap 
will exceed y tends to 1 as r — > oo. Thus finite configurations of more than 
one particle can carry only zero probability in any quasi-stationary measure. 
Of course, a measure with exactly one particle is quasi-stationary. □ 

5.2. Some auxiliary estimates. Given an initial configuration oj = {x n }, 
the probability distribution of the position of the leading particle after r > 
steps is dP^\x), with 

Pw\ x ) = Prob({at time r all particles are on (— oo,x]}) 

(5.8) 

= Y[[l-P(S T >x-x n )}. 

n 

We shall need to compare dP& (x) with the probability distribution associ- 
ated with the function 



(5.9) PM (x) = expj - £ P(S T > x - x n ) j = e 



-F u;t (x) 



Remark. It is instructive to note that dP^ T \x) is the probability distri- 
bution of the maximum of a modified process, in which at first each particle is 
replaced by a random number of descendents, with the Poisson distribution 
p n = e _1 /n!, and then each particle evolves by r independent increments, as 
in the Indy-500 model. Conditioned on the starting configuration, the mod- 
ified process is (instantaneously) a Poisson process. The probability that its 
maximum is in (— oo,x] is given by 



e- 1 



n 

n 

(5.10) 

= exp 



^ (l-P{S T >x-x n )) n 

^ n! 

n 

J2P(Sr>X-X n ))=P^(x). 
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Our first goal is to show that the probability measures dP^ T \x) and 
dPu {x) are "typically" — in a suitable stochastic sense — asymptotic to each 
other as r — > oo. This statement is not true for some uj, and it is not difficult 
to construct examples of configurations for which it does not hold. We note 
that it is easy to show that the step described by the graph of P^ (•) remains 

tight, in the sense that the width of the intervals {x:5 < P^(x) < 1 — 5} 
does not spread indefinitely, as r — ► oo. 

Lemma 5.3. For any quasi-stationary measure 



(5.11) eJ r supP(S T >x-x n )dPW(x) 

\J—oo n 

Furthermore, 

(5.12) eJsu P \P^(x)-P^(x)\] — 1 0. 



— ► 0. 

r— >oo 



r^oo 



Remark. The supremum in (5.11) is clearly attained at n = 1 (by mono- 

tonicity). Since dP£\x) is a probability measure, and the c.d.f. of S T is a 
bounded function, the statement means that the maximum typically occurs 
in a region whose a priori probability of being reached by any specific point 
is asymptotically zero. 

Proof of Lemma 5.3. Due to the spreading property of convolutions 
of probability measures (see [4], Lemma 9.4.1), for any D < oo 

b(r,D) = sup P(x<S T <x + D) — ► 0. 



T — >0O 



Observe that P^\x) < P^\x) < 1 for all x. Let us pick A > such that 

e _x(l+Ax) < 1 _ x Vx€[0,i]. 

Thus if P(S T > x) < we have 

(5.13) PW(X) < p(T)( X )l/lW(Sr>x)]. 

Suppose that x is such that P(S T > x) <e. Then 

Pu T \x) -P^\x) < sup \u 1/(1+Xe) -u\. 
«e[o,i] 

Suppose that x is such that P(S T > x — x±) > e. 

Let no = - In -. Then for all r large enough and for all n < no, 

b[r,-x n ) < -. 
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Consequently, 

P(S T > x - x n ) > P(S T > x) - 6(r, x n ) > e - . 

Then 

(2/ e )ln(l/e) 

-Y. P (Sr>X-X n )<- J2 n<" ln - 
n n =0 6 

and therefore 

(5.14) P>)(x) < e - ln (i/^) < £ 
So in this case we obtain 

PW (*)<£. □ 

Putting the above together, we have: 

LEMMA 5.4. If be a quasi- stationary measure, then for each e > 0, 

(5.15) M ({^:dist(dpW,dpM) > e}) — 0, 

r — >oo 

where dist is the distance between the two measures, defined as 

(5.16) dist(dP, dP) = sup j J h(x) dP{x) - J h(x) dP(x) 

Proof. The distributions dPu{x) and dPu T \x) can be written as 
(5.17) 

rfPi T) (*) = E ^ X " X k) X eX p| " E P ( S r >X-X n )\. 

k { n J 

By Lemma 5.3 we obtain that 

\P^ T) (x)-P^ T) (x)\<e Vx. 

If P(<SV > x) < e, then we obtain by the same arguments as in the previous 
lemma that 

\dP^{x)-dP^{x)\<edP^{x). 
Integrating with respect to tJt\ — over the x such that P(S T > x) < e, we 

|| " 1 1 oo 

obtain that the result is small. 
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If P{S T >x)>e, then 

Y[[l- P(S T >x- x n )]<e and expi - ^ P(S T > x - x n ) I < e. 
Consequently for such a? 

\\h\\co u { ) 

<V~e J^dP{S T >x-x n )e^-]^ P ( S r>x-x n )^ 

< const y/e, 

and also using 1 — x < e~ x for x > 

< const y/e. □ 

5.3. "Poissonization" — the proof. We are now ready to prove the main 
result of this section. 

Proof OF Theorem 5.1. Due to the quasi-stationarity of the measure 
H, one may evaluate G^{f) by taking the average of the future expectation 
value of exp{— £)n/(fi — Vn)}i corresponding to the configuration uj as it 
appears at time t = 0. 

In the following argument we fix the (nonnegative) "test function" /, 
and take D < oo such that supp/ C [— D,0]. In the approximations which 
follow we use the fact that exp{— J2 f(yi — Un)} is a bounded function (< 1), 
which is integrated against a probability measure. As before, w-dependent 
quantities are denoted o(l) if in the limit r — > they tend to "in law," 
that is, the probability distribution which they inherited from u is nonzero 
only for [0, e] for any e > 0. 

The conditional expectation of the future value of exp{— J2f(lJi ~ Un)}, 
conditioned on the initial configuration u, is 



poo 

(5.18) = / e- f{0) dP(S T >x-x k ) 
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X [] [1 - P(S T > X - X n )] 
n^k 

too^ S{X ~ V) dP{S T >y-X n ) 



X 



n 



(l-P(S T >x-x n )) 

where dP{S T > x — x^) is the probability that the fcth particle is at x at 
time r, rin^fcll ~~ P(S T > x — x n )] is the probability that other particles 

r e-f( x -») dP(S T >y-x n ) 
are at (— oo,x] at time r, and — °° nzrprs > x -x )) ^ s the expectation of 

e~f( x ~ yn > given that the particle which is at x n at t = is at (— oo,x] at 
time r. 

As in the previous discussion, the contribution of x such that P(S T >x)> 
e to the integral in (5.18) is negligible. 

Consider x such that P(S T > x) < e. We can write 

jj f^e-te-yUPjSr^y-Xn) 
11 l-P(S T >x-x n ) 



(5.19) = J] 

nj^k L 



1 _£ sp (l-e-n*-v))dP{S T >y-x n ) 



1 — P(S T > x — x n ) 



(l + o(l))exp{- f (l-e- f ^)d(Y,P(Sr> 



As noted in (5.3), the normalizing shift has no effect on Gf(J)- The result 
is the first of the two equations in (5.6). The second equation is an immediate 
corollary of the first one, since 

(5.20) g*F u . T = F u . T+1 . □ 

For a later use, let us note that the arguments used in the above discussion 
readily imply the following two bounds. 

Corollary 5.5. For any e > 0, there is W(e) < oo such that 

(5.21) eJ f d[e-^ Fu ^]]<e 

\J\x\>W(e) J 

and 

(5.22) eJ f d[e- g ** rF »^ x) ]) < e. 

\J\x\>W(e) ) 

Proof. Let / = I[o ) w(e)]- Denote 

(5.23) <l>(W(e)) =E fl [e- 1 ^ w ^)M-y^]. 



16 A. RUZMAIKINA AND M. AIZENMAN 

Since I[o,w(e)] ( x ) 1 for i G 1 and since, in a typical configuration, 

the number of particles within distance W{e) behind the leader increases 
to oo as W{e) increases, <fi(W(e)) must decay monotonically to as W(e) 
increases. By taking / = I[o,w(e)]> we see that 

4>(W{e)) = I Li{duj) r de- MF <"^ 

J J —oo 

(5.24) 

x e -(l-e- 1 )(ATF W!T (x-W'( e ))-JVF (1 , iT (*)) + 0(e T ). 

We can get an estimate on AfF w - T (W(e)) from (5.24) by restricting the 
range of integration from W(e) to oo and using that J\fF u . T (x — W{e)) — 
MF UJ]T (x) < 1. Then, for x > W(e) we obtain 

<KW(e)) > I n(du) r de~ NF ^ x) 

J JW(e) 

(5.25) X e Hl-e-^){NF^Ax-W{e))-NF^ T {x)) + 0(e t ) 

> e-^ 1 " 6 " 1 ) ffi(du)(l - e-^W)) + (s T ). 



Similarly, by restricting the range of integration from to oo and using 
that NF^x - W{e)) - NF w . T {x) < N F U]T {-W {e)) for x > 0, we obtain 

<f>(W(e))> J Kdu) J™ tfe-^frC*) 

(5.26) X e -(^e-^F u . T (x-W(e))-^F u;T (x)) + 

> (1 - e" 1 ) / ^du)e^ l ~ e ^ NF ^- w ^ + 0(e T ). 



Equations (5.25) and (5.26) prove the first part of the corollary. 

To prove (5.22) we observe that from the previous part it follows that for 
all t large enough, and sufficiently large W(e), 

E„ / de- NF ^+Hx)<-. 

J\x\>W{e)/2 2 

Since for sufficiently large W{e) and u in a set of measure 1 — | , 

z < W{£) 

we obtain that 



J\x\>W(e) 
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□ 

Corollary 5.5 will be used for an approximation of Gj^fp (f) by a quantity 
which has better continuity properties as a functional of F. 

6. The Poisson density as a Laplace transform of a random positive 
measure. We shall next show that the quasi-stationary measure p can be 
presented as equivalent to a random Poisson process whose density is the 
Laplace transform of a random positive measure on K. [Due to the invari- 
ance of B under uniform shifts, with no additional restriction the measures 
may be adjusted so that p(R) = 1.] 

Let A4 be the space of finite measures on [0, oo). To each p € M. we 
associate the Laplace transform function 

(•OO 

(6.1) R P {x)= / e- xu p(du). 

Jo 

We denote by Tl the space of such functions, that is, JF L = {R p (-)\p £ Ai}. 

We shall need to consider "ensemble averages" over randomly chosen el- 
ements of J\A. These are described by probability measures on J\4, which 
would always be understood to be defined on the natural cr-algebra on J\4, 
for which the measures of intervals, p([a,b]), are measurable functions of p. 
Our goal in this section is to prove the following statement. 

Theorem 6.1. Under the assumptions of Theorem 4.2, there exists a 
probability measure, v(dp), on M such that for any compactly supported 
positive function f on M, 

(6-2) G„(f)= f u(dp)G Rp (f), 

Jm 

and furthermore, 

(6-3) G,(f)= f v(dp)G Rp * 9 (f). 

J M 

For Laplace transform functions F = R p , shifts correspond to transforma- 
tions of the form 

(6.4) p{du) ^ e - au p(du), 

and the normalization condition (5.1) corresponds to p(M) = 1, that is, 
p G M being a probability measure. In view of the invariance (5.3), this 
normalization condition may be freely added as a restriction of the support 
of v(dp) in the statement of Theorem 6.1. 

While the result presented in the previous section required only quasi- 
stationarity, we shall now make use of the additional assumptions listed in 
the main theorem (Theorem 4.2). 
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In the derivation of Theorem 6.1 we shall apply what may be regarded as 
the principle of the equivalence of ensembles, in the language of statistical 
mechanics. Specifically, we need the following result, which, as is explained 
in the Appendix, is a refinement of the "Bahadur-Rao theorem" of large 
deviation theory. 

Theorem A.l. Let u\,U2, ■ ■ ■ be i.i.d. random variables with expectation 
E g u and a common probability distribution g(u), which has a density and a 
finite moment generating function, J e iru g{v) du = < oo for all n. Then, 
for any < K < A'(oo) and < 8 < i there is £ t k b — > such that for all 

q 6 [E g u, K] and \x\ < , 

, , Prob({ni + u 2 + • • • + u T > x + qr}) _ 
Prob({wi + «2 H Vu T > qr}) 

with rj = rj(q) determined by the condition 

f ue vu g(u) du 



(6.6) 



feWg{y)dy 



In our analysis we shall need a bound on the front velocity, and on the 
possible propagation of particles from the far tail. 

Lemma 6.2. Let fi be a quasi-stationary g-regular measure with a density 
satisfying the assumptions (4.2) and (4.3) of Theorem 4.2. Then: 

(i) For any r large enough, for to in a set of measure 1 — e, 

S f 

(6.7) z w - T < —t + const where S = In / e 2Xx g(x) dx. 

2A J 

(ii) There exist a M (M) and 8 9 {t) such that the probability of the comple- 
ment of the event 

A t -d,k,m = : the configuration obtained after r steps will have not more 
than M particles with y n > y\ — D, and all of them made 
a total jump less than Kt + z Wjr — x n in time from to r} 

satisfies 

(6.8) Viob(A c T . DtKM ) < a^M, D) + 8 g {r) + C g>fl e-^ K - K °> 
with a„(M,D) — > for each D < oo, B„(t) — > and 5 > 0. 

Remark. In the proof below we shall apply the last bound in the double 
limit: lim^^oo lim r _ KX) , with M chosen so that 1 <C M <C K. 
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PROOF of Lemma 6.2. (i) By (4.3) and Markov inequality, 

^({t)(-x„)< ? n}>e 2Am )< e - Am . 

Therefore by the Borel-Cantelli lemma, 

P^m-xn) < m} > e 2Xm i.o.} = 0. 

This implies that for any e there exists niQ such that on a set of u of measure 
1 — e, {#(— x n ) <m}< e 2Xm for all m > m . 
Using the definition of F U] q we obtain 

(6.9) F^ {x) < e -2Amm(x,- mo ) Vx<0. 

Therefore, 
(6.10) 

< const J e~ 2X{x - y) g { * T) (y) dy < const e ~ 2Xx+TS . 

For x = |j + const we thus obtain F LU - T (x) < 1. It follows by definition that 
Zw,t < fx + const. 

(ii) The probability that the first condition does not hold in the definition 
of A t -d : k,m is, by the quasi-stationarity of fi, 

an(M,D) = n(u> :more than M particles 

(6.11) 

are within distance D of the leader at t = 0). 

This quantity vanishes for M — ► oo because the number of particles in [yi — 
D, y\] is almost surely finite. 

To estimate the remaining probability of the complement of the event 
A T] d,k,m we split it into two cases, based on the distance which the front 
advances in time r. That distance is at least the total displacement of the 
particle which is initially at 0. The probability that this displacement is less 
than (E g u — l)r is dominated by the quantity 

Prob(^, r < {E g u - l)r) < Prob(S T < (E g u - l)r) = p g (r). 

The choice of 1 is somewhat arbitrary, but even so, standard large devia- 
tion arguments which are applicable under the assumption (4.2) imply that 
j3 g {r) decays exponentially. 

The contribution of the other case is bounded by the probability of the 
following event: 

Prob(at least one of the particles of u will advance in r steps 
a distance greater than [— x n + (E g u — 1 + K)t\) 
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< E M Prob(S T > -x n + (E g u - 1 + K)r)j 

< (^E g {e aST )e- a[ - Xn+iE ° u - 1+K)T] ^J 

< J e a{u - E v u) g{u)due- a{K - l) ^ J ^V^-^J , 

where a > is an adjustable constant. The last factor is finite for < a < A 
since under the assumed exponential bound (4.3), 



(6.13) 



^E^"^) = E fl (a J dye- a ^l[y>-x n yj 
<aj dye- ay Ae Xy 



a — A 



The claimed estimate readily follows (choosing A < a, and defining 5 > 
correspondingly). □ 

Proof of Theorem 6.1. Applying Theorem A.l to the function de- 
fined by (5.4), we find that 



NF^ T {x) = £P(SV > X + Z^r - x n ) 

n 

= y ^ P(S T > z u ^ T — x r 

-K(e)r<x n <0 



P(S T ^ x -\- Ztjj-^j- x Ti 
P(S T ^ z w ^ T x n ) 



(6.14) 



+ £ P ( S r > Z w,T ~X n + x) 
x„<—K(e)t 

£ P ( S r > z»;t ~ x n )e-^ z ^-^l^ x [l + 0{e T )\ 

-K(e)r<x n <0 



^ %UJ,T X n + 

Xn<—K{e)r 

= p^ T {du)e- ux [l + 0(e T )} 
Jo 

x„<— K{e)r 

with Puj- T (du) defined as the probability measure with weights P(S T > z w;T 
x n ) at the points i]{ ZuJ ' T Xn )- 
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We will now estimate the remainder term X^x„<— k(s)t P(^t — z w,t ~ x n + 
x). 

In the case when lim^—nx, A' (77) < 00 (in the case when the supremum of 
the support of g(x) is finite), the remainder term is zero for large K(e) [e.g., 
if K(e) > A'(oo) and x = 0(t p )). 

In the case when lim^oo A' (77) = 00, the remainder term can be estimated 
using the large deviation arguments. By using (A. 2) in Theorem A.l and 
(6.9) we obtain 

y ] P(S T ^ %u),t — x n + x) 
x„<—K(e)t 



POO 

< / P{m + ---+u T >y + z 0J , T + x)e 2Xy dy 

JK(s)t 

(6.15) = / expj-rA* 

JK{e)r I V T 



jT ex P { („ {y±^±^ ) t } dQ ^M) {t) 



e 2Xy dy 



= 0(e T ). 

The last equality in (6.15) follows because by convexity of A* 

A I ■ J» for all y > K (e)r, 

and because the factor in the square brackets in (6.15) is small. 
Therefore 

POO 

(6.16) NF u) , T (x)= e~ ux PuJ . T (du)(l + 0(e T )). 

Jo 

We observe that 

\tfF urr (x)-R p {x)\<eR p (x), 
\AfF^ T (x)-R' p (x)\<eR' p (x), 

(6.17) 

\g*MF ul . T (x) -g*R p (x)\ <eg*R p (x), 
I (9 * MF u . r )'{x) - (g * R p )'(x)\ < e(g * R p )'(x). 
Using (6.17) and Corollary 5.5 we obtain 

/W(e) ( rx *\ 

de- R »Mex V \- (l- e -/(*-w))(-dR ( y ))L+ e 
-W(e) I J-oo J 

(6.18) =G wie)jRp (f) + s 



22 



A. RUZMAIKINA AND M. AIZENMAN 



/W(e) ( rx > 

d e -9* R ptoexp\- / (1 - e- f( - x - y ^)(-dg * RJy))\ + e 
-W(s) I J-oo ) 



rW(e) 
Gw{e),g*R p (f) + e 



[Equation (6.18) will serve as a definition of Gw(e),R p (f)-] 

From (6.16) we observe that for all u in a set of measure 1 — e and for 
every K S> 1 there exists a i^i 3> 1 depending on i"T such that 

poo 

(6.19) / e Du Pu , T {du) < V P(S T >x + ^ ;r -x n ). 

h(E a u+K) Xn <- Kl r 

We can choose i^i by requiring that for all x n < —Kit, 

! < -Eg« + K, 

T 

for example, K\ = E g u + K — and we used (6.7). 

From Lemma 6.2 [see (6.8)], applied with M = \f~K (or any other choice 
with 1<M< K), we find that under the assumptions listed above, for any 
D < oo there exist ed(K) with which 



limsupE^I / e Du p UJ - T (du) 

t^oo \Jr)(EgU+K) 



(6.20) 



<limsupEj P(S T >z u . T -x n -D))+s D (K) 

T ^°° \x n <-K lT / 



e D (K)-^ 0. 



The correspondence lo \— > p^ T defines a mapping from the space of con- 
figurations £1 into the space of measures on R, with values restricted to 
the subset of probability measures. Corresponding to this map is one which 
takes the measure /ionS] into a probability measure on A4 which we shall 
denote by v T . By this definition, for any measurable function X : M — > R, 



(6.21) / X(p)v T (dp) = / X(p w . iT )v(duj). 

The space of probability measures on compact subsets of R is compact, 
and so is the space of probability measures on this space. While we do 
not have such compactness (since the measures of A4 are defined over the 
noncompact R), (6.20) with any fixed D > implies that the sequence of 
measures v T is tight and that it has a subsequence u Tn which converges in 
the corresponding "weak topology" as r n — > oo. Let v be a limit of such a 
subsequence. [To prove the tightness of u Tn we observe that it is possible to 
show that for all r, R Pul . T (x) < M(x) for some function M(x) except for lo 
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in a set of measure e. The set of p for which R„(x) < M(x) is compact.] We 
claim that for every positive / of compact support, supp/ C [— D,0], 



M 



[G W (e),R p (f)Hdp) + e 



(6.22) 



hm / [G wie)tRp (f)]u Tn (dp)+e = G ll (f), 



M 



[Gw(e),g*R p {f)]v{dp)+£ 



Jim / [Gw(e), g *R p {f)]vT n (dp) + e = G M (/). 



The weak convergence means that for any continuous function <Y : .M 



*(p)K^) = hm / X(p)v Tn (dp) 
(6.23) ^ n °° Jm 



The continuity argument does not apply immediately to the function which 
we are interested in: 



Gw(e),R p (f) 

/W{e) ( rx 

d[e- R »M]e W \- [l-e-H x -ri]d(-R p (y))\, 
-W(e) I Jx-D J 



^ We) 



-W{e) 

which is not continuous in p. However, G-w( £ ),r (/) can be approximated 
arbitrarily well, in the appropriate L\ sense, by functionals which are con- 
tinuous. 

The function Gyju\ R (/) is not continuous in p. The difficulty is that 
R p (x) can be affected by small changes in the measure p if those occur at 
high values of the Laplace variable u. However, we do obtain a continuous 
function by replacing R p in (6.24) by R\ KP with 

(6.25) iKp(du) = l [0M E a u+K)} 0) • p(du) . 

It is easy to see that 



e- xu p(du) )v{dp) 

M \Jr){E 3 u+K) J 

(6.26) <limsup / ( f°° e~ xu p^du) J fi(du 

r^oo Jn\J v (E„u+K) ' ) 



<e x (K)-^ 0, 

K— >oo 
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where the first inequality is by the generalized version of Fatou's lemma, 
and the second is by (6.20). 

Due to the fact that / is compactly supported and the integration in x is 
over [— W(e), W(e)], the difference 

\G W (e),R p (f) ~ G W{£ ), R (f)\ dv{p) 

M 

is affected only by values of x G [— W(e) — D, W(e)]. 

Taking x in this interval, we observe that (6.26) implies that 

roc 

/ e- ux p(du) < e, 

Jr](E g u+K) 

except on the set of u of measure e. 

The difference J M \G w ^ jRp (f) - G w ^ >RlKp (f)\ dv[p) is controlled by 
\Rp{x) — R\ KP {x)\ and by \R' p (x) — R[ K p (x) \ , which are small for x G [— W(e) — 
D,W(e)] except on the set of p of measure e, since J^E g u+K) e ~ UX p(du) is 
small. 

One can verify by standard arguments that for K finite Gw( £ ),Rx p (f) 
is continuous in p, and that this continuity and the approximation bounds 
listed above imply (6.22), thereby proving the first part of Theorem 6.1. The 
second part is proved via similar arguments. □ 

7. Monotonicity arguments. In this section we develop some monotonic- 
ity tools, which will be applied to prove that if a measure p, has the properties 
listed in Theorem 6.1, then the corresponding measures p are z^-almost surely 
concentrated on points, that is, the Poisson densities R p are almost surely 
pure exponential. 



7.1. The contraction property of convolutions within T^. The space 
whose elements are positive decreasing continuous functions on M, is par- 
tially ordered by the following relation. 

Definition 7.1. For F, G G T we say that G is steeper than F if the 
level intervals of G are shorter than those of F, in the sense that for any 
< a < b < oo, 

(7.1) (0 <)G'\a) - G'^b) < F-\a) - F~ 1 (6). 

We adapt the convention that for the (monotone) functions G G T the in- 
verse is defined (for a > 0) by 



(7.2) 



G (a) = inf {x G M : G(x) < a}. 
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It is easy to see that, within the class of monotone functions J-, an equiv- 
alent formulation of the relation "G? is steeper than F v is that for any u > 0, 



(7.3) 



G(x)=F(y) 



G(x + u) <F(y + u). 



Also equivalent is such a principle with the reversed inequality and u < 0. 

Of particular interest for us is the subspace J~l of Laplace transforms of 
positive measures. We shall show that within this space, the convolution 
with a probability measure g(x) dx makes a function steeper. (It is shown 
below that the appropriately shifted R p * g is in J-l-) A key step towards 
this result, which is also of independent interest, is the following lemma. 

Lemma 7.2. Let F = R p F £ J-l satisfy the normalization condition F(0) = 
1 (i.e., F = J\fF), and let G = R p c £ Tl be related to it by 

(7.4) G = M(F * g) 

for some probability measure g(x) dx. Then, for all A > 0, 



(7.5) 



A r\ 

p G {du) < / p F (du). 
o Jo 



Proof. The relation between F and G is such that for some normalizing 
constant z£l, 



/o 
-c 

with £(•) defined by 

(7.7) 

Thus 
(7.8) 



-(x—y+z)u F 



/ (du) 



g(y) dy 



-xu S(u) F 



p r {du) 



D S(u) 



p G {du) 



g{y)dy. 



,S(u) F 



p (du), 



and the normalization conditions F(0) = G(0) = 1 imply 



(7.9) 



S{u) n F 



p b (du) 



p F (du) 



10 JO 

The function S(-) is convex, which is easily verified by showing that S" > 
0, by general properties of integrals of the form (7.7), and satisfies 5(0) = 
(since g(x)dx is a probability measure). It has, therefore, to be the case 
that either p F (du) is concentrated at a point (where S = 0), or else S(-) 
is negative on [0,u) and positive on (u,oo) for some u > 0. The claimed 
concentration statement (7.5) is obviously true for all A S [0,u]. For A > u, 
we note that 



(7.10) 



,S(u) F 



p P (du)> I p p (du). 
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By subtracting (7.9) from (7.10), we find that the claimed (7.5) is valid 
also for A > u. □ 

Theorem 7.3. For any F = R p G Tl and a probability measure g(x) dx 
on R ; the function M(F * g) is steeper than F . 

Proof. Our goal is to derive the inequality (7.1) for G = M(F *g) (and 
a <b). By simple approximation arguments, it suffices to do that assuming 
linx^.ooF^) = oo. 

We claim that 

Af(F*g)(x) <F(x) forx>0, 

(7.11) 

Af(F*g)(x) > F(x) forx<0. 

We find that the functions F and G = N(F * g) are related just as in the 
previous lemma. In order to convert the concentration statement (7.5) into 
one relating G(-) with F(-), we write, using Fubini's lemma (or integration 
by parts), 

roc r pX 

forx>0 M(F*g)(x)=x d\e~ Xx / e s{u) p{du) , 

Jo Uo 

roc r /"oo 

(7.12) forx<0 N(F*g)(x)=x d\e~ Xx / e s{u) p{du 

Jo Ux 



3C 



+ / e s Mp(du), 
Jo 

with the corresponding relations holding for F without the factors e s ( u \ 
The inequalities (7.11) follow now by inserting here the relations (7.5), 
(7.10) and (7.9). 

We note that if F and N(F * g) were shifted so as to be equal at a different 
value of x, then the argument above would also go through. Therefore we 
obtain that M(F * g) is steeper than F. □ 

The partial order "G is steeper than F n is preserved when any of the 
functions is modified by a uniform shift, and also when each is replaced by a 
common monotone function of itself, for example, {F, G} replaced by {1 — 
e~ F , 1 — e~ G }. Following is a useful property of this partial order (another 
one is presented in Appendix A. 2). 

Lemma 7.4. Let F,G€ T be continuous and strictly monotone decreas- 
ing functions with lirn^-oo F(x) = \im x ^ > - 00 G(x) = oo. If G is steeper than 
F, then, for any u > 0, 

(7.13) f e -[G(x-u)-G(x)] de -G( X ) < f e -[F(x-u)-F(x)] de ~F(x) _ 
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Furthermore, the inequality is strict unless G is a translate of F (and vice 
versa). 

Proof. The statement is a simple consequence of the following formula, 
and (7.3): 

l e -[F(x-u)-F{x)] de -F(x) _ I e -[ G (x-u)-G(x)] de ~G(x) 

(7.14) 

dz r e - J!, ( J!, ~ 1 W-") - e-GiG-^zyuh 



□ 

An additional result related to this notion, which may be of independent 
interest, is presented in Appendix A. 2. 

8. Proof of the main result. We shall now apply the monotonicity argu- 
ments for the last leg of the proof of our main result. Theorem 4.2 is clearly 
implied by the following statement (see Theorem 6.1). 

Theorem 8.1. Let p be a measure on the space of configurations J7 ; 
which admits a representation as a random Poisson process, described by a 
probability measure v(dp) on M. as in Theorem 6.1, for which both (6.2) and 
(6.3) hold. Then the support of the Laplace measure dp is v-almost surely a 
point; that is, the functions R p are almost surely pure exponentials. 

Proof. Let us consider the probability that the first gap exceeds some 
u > 0. For a Poisson process, a simple calculation yields 



. . \ / oo 

(8.1) jg (Poisson) ^ _ X2 > u)= e - F ^- u \-dF(x)). 

J — oo 

Therefore, 

E M (aci - x 2 > u) = f n(du>) (°° e- F< - x - u) (-dF(x)). 

Substituting this in (6.2), or in (6.3), one obtains the corresponding expec- 
tation for the measure ji. Subtracting the two expressions, we find that 

p |~ roc roc 

(8.2) 0=/ v{dp) / e- R "^- u UR p {x) - e -Rp*g(*-«) dR p * g(x) . 

By the analysis in the previous section (Theorem 7.3 and Lemma 7.4), the 
difference in the square brackets in (8.2) is nonnegative. Thus, this relation 
implies that 

f e -R P (x-u) dRp ^ _ J e -R p *g{x-u) dRp ^ g ^ = 

(8.3) 

for z/-almost every p. 
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Furthermore, by Lemma 7.4 the equality yields that i^-almost surely R p * g 
coincides with one of the translates of R p . The only functions {F = R p ) with 
this property in Tl (or for that matter in T\ see [3]) are pure exponentials, 
which correspond to p concentrated at a point. □ 

APPENDIX 

A.l. Useful statements from the theory of large deviations. Our goal 
here is to derive Theorem A.l which was used in Section 6. Its statement may 
be read as an expression of the "equivalence of ensembles" — in statistical 
mechanical terms. The following notation will be used in the theorem. 

A(A) = lnE[e Aui ], A*(y) = sup(Ay - A(A)). 

A 

The result we used in Section 6 is: 



Theorem A.l. Let u\,u%... be i.i.d. random variables with a common 
probability distribution g, which has a density and a finite everywhere mo- 
ment generating function. Then, for any < K < A'(oo) and < (3 < 1/2 
there is e t -k 8 — ► such that for all q G \E„u,K] and \x\ < t$ , 

v Prob({m + u 2 + • • • + u T > x + gr}) _ 
Prob({iii + u 2 H h u T > qr}) 

with rj = 77(g) determined by the condition 77(g) = A*'(g). 

Proof. We will assume that E g u = 0, since we can replace the random 
variables Ui by — E g u. We will use the same notation as in the proof of 
the Bahadur-Rao theorem (see [5]). We denote 



« r U|2)) = ,(2) ) /rA.(,(2 

Ui-y/r 



Yi 



. Y X + ---+Y T 



and consider a new measure P defined by its Radon-Nikodym derivative 

aip(v(y/-r)) 



dP 



■{x) = e 
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Let also denote the distribution function of W T with respect to 

pWn/T)) _ j g eaS y i show then that Y{ are i.i.d. with mean and variance 
1 with respect to p(w/ T )) . Therefore W T has mean and variance 1 with 
respect to Q^ (s//r)) . 

By analogy with the proof of the Bahadur-Rao theorem (see [5] ) , we can 
write 

r OO 



(A.2) P(u± + --- + u T >y) = e 



-tA*(v/t) 



-i'riviv /r))t (IqMv/t)) 



For further consideration, we need to estimate the ratio 
P(u\ + u 2 H h n T > x + y) 



(A.3) 



P(ui + U2 H h^T > y) 



-tA* ((i+s)/t)+tA* (j//t) 



t) Jo 



00 P -^r(»7((a!+y)/T))t jq^ x +v)I t )) 



J™e-Mv(yM)t dQ ( T v(yM) (t) 
By using Taylor's expansion we can estimate the exponent in (A.3): 



(A.4) 



A* 



x + y 



+a*( y~ 



-1-20 



where, to estimate the remainder term in (A.4), we use the integral form 
of the remainder in Taylor series and that - < K, \x\ < t@ , A*" = 7^77 < 00, 
convexity of A and the assumption that the Laplace transform of g is finite. 
It remains to show that the prefactor in (A.3) is 

r°° e -Mv((x+v)H)t dQ Mix+y)/r)) (t) 
(A. 5) r(x,y) = ^ . , . — = l + 0(e T ). 

By the Berry-Esseen theorem (see [7]), 



sup 

X 



-t 2 /2 



dt 



< 



33 Euf 1 



Therefore, 



(A.6) 



4 (Varui) 3 /2 

ip T M(x+y)/T))t dQ^ x+y ^ T ^(t) 



o 



-<p T (ri((x+y)M)t-t 2 /2 



"7? 



2tt 



■dt 



x + y 



0(V 



This formula is especially useful when ip T < O(l) (i.e., when 77 is small) and 
the first term on the right-hand side of (A.6) is much larger than the second 
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term. In this case we obtain 

_ S?e-*TM*+v)M)*-*/*dt + 0(l/SF) 
{ ' { ,V) f™e-^(r>(y/T))t-ty2 dt + 0{1/ ^ T) ■ 

If y is such that 0(1) <ip T < (^(t 1 / 2 ), we write the integral as 

oo ;>oo 

e -Mv(vM)t d Q(v(y/r)) = / e -Mn(v/r))t q {p(y/r)) ^ dt> 

J 

where g T is the density of Q T . By the analog of the Berry-Esseen theorem 
for densities (see [7]), 



'( — ^ ) 



(A.8) sup 

From (A.8) we obtain 

r(x v)= [°° e -Mri((x+v)H)t-ty2 dt + 1 o(— 

(A.9) _ x 

x ( r e -^(v(yM)t-t*/2 dt+ [ (±: 

uo Vr{v{y/T)) \V T 

The proof of (A. 5) now consists of showing that 

e -^ T (r ? ((x+y)/r))^t 2 /2 df _ f°° g -Vv (v(yM)t-t 2 /2 ^ 

io 

''OO ~ 

e -i>r(v(yM)t-t 2 /2 dt 



fC 

< 0(0 / 

JO 



□ 



A. 2. A class of monotone functionals over Since the notion intro- 

duced in Section 7 may be of independent interest, let us present here a 
related result, which offers another instructive insight on the contraction 
properties of convolutions in T . 

Theorem A. 2. Let F,G G T with G steeper than F. Then, for any 
positive and continuous function *S> on [0, oo) which vanishes at and oo, 

/OO /"OO 
dm(G{t)) < / dtV(F(t)). 
-oo J — oo 

Furthermore, if ^ is strictly positive on (0, oo), and G and F are both left- 
continuous, then the inequality is strict unless G is a translate of F. 

PROOF. By standard approximation arguments (e.g., using local approx- 
imations by polynomials), it suffices to establish (A.10) under the assump- 
tion that ^ is piecewise strictly monotone. 
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Employing Fubini's lemma, or Lebesgue's "layered cake" formula for the 
integral, 

/oo poo poo 

dtV(F(t))= dX dtI[$(F(t)) > A], 
-oo JO J —oo 

Under the added assumption on ^, the set {f Gl: ^(F(t)) > A} is a union 
of level intervals of F, of the form {t G R:a/(A) < F(t) < bj(X)} (with 
{[dj(X),bj(X)]}j determined as the level sets of > A}). 

The integral over t on the right-hand side of (A.ll) produces the sum 
of the lengths of the level-intervals of F. When F is replaced by G, the 
corresponding intervals can only get shorter, since G is assumed to be steeper 
than F, and thus (A. 10) holds. 

In view of the above, the conditions for the strict monotonicity sound 
reasonable. However, since the strict monotonicity is very significant it may 
be instructive to make the argument explicit. (What follows makes the ar- 
gument given just above redundant; however, we keep it because of its sim- 
plicity.) It is convenient to rearrange the above argument as follows. Using 
our convention for the inverse function, 

/oo poo poo 

dm(F(t))= dt{F)V{F)= dF~ 1 (a)^(a), 
-oo J — oo J — oo 



and thus 



(A.13) 



/oo 
dtV(G(t)) 
-oo 

/oo 
[dF^ia) -rfG" 1 (a)]^(a), 
-oo 



where dF~ 1 (a) and dG~ 1 (a) are measures on M. The assumed relation (7.1) 
implies that the difference dF~ 1 (a) — dG~ 1 (a) is itself a positive measure. 
The vanishing of its integral against ^ is therefore possible only if this 
measure is supported in the set ^ _1 (0), but that set (viewed as the set of 
values of the functions F and G) contains at most the boundary point a = 0. 
It follows that if the inequality (A. 10) is saturated, then the two Stieltjes 
measures are equal in (0, oo), and thus 

(A.14) ^ _1 (a) -G" 1 (a) = Const, 

which means that F and G differ by a shift. □ 



This implies another monotonicity principle, which expresses the fact that 
convolutions make functions in Tl steeper. 
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Corollary A.l. For any function F £ Tl and a probability measure 
g(x) dx on R, 

(A. 15) E^ lsson) (x n - x n+ i) < E^ >olsson) (x n - x n+1 ) for all n > 1, 

and the inequality is strict unless either both quantities are infinite, or F(x) = 
e -s{x-z) j Qr some s > o and z £ R. 

Proof. The mean value of the gap may be computed with the help of 
the expression 

/oo 
{I[t > x n+ i) - I[t > x n ]}dt. 
-oo 

A simple calculation yields 

(A.17) E^ oisson \x n - x n+l ) = r dW n (F(t)) 

J —oo 

with fy n (F) = e~ F ^ . Theorem A. 2 applies to such quantities. □ 

We did not base the proof of Theorem 6.1 on this observation [i.e., use 
in Section 7 (A. 15) instead of (8.1)] since this argument is conclusive only 
when the above expected value is known to be finite for some n < oo, and 
we preferred not to limit the proof by such an assumption (and had no need 
to). 
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